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ABSTRACT 

The most detailed constraints on the accelerating expansion of the universe and details of 
nature of dark energy are derived from the high redshift supernova data, assuming that the 
errors in the measurements are Gaussian in nature. There is a possibility that there are direction 
dependent systematics in the data, either due to uncorrected, known physical processes or because 
there are tiny departures from the cosmological principle, making the universe slightly anisotropic. 
To investigate this possibility we introduc e a sta tistic based on the extreme value theory and apply 
it to the gold data set from iRiess et al. ( 2004 ). Our analysis indicates a systematic, direction 
dependent non-gaussianity at about one sigma level. 



Subject headings: cosmology: cosmological parameters 
— supernovae: general 

Introduction 



During the last decade the possibility that the 
expansion of our Universe is accelerating has been 
put on a firm footing. The c ombined analysis 
of high redshift supernova dat a dPerlmutter et aT 



19991 : IRiess et al.l Il998l 120021 12004) . along with 
observ ations of cosmic microwave background 



(jBenoit et al.ll2002l;lPage et al.ll2006l) or large scale 



structure (Pe rcival et alj 2002; iTegmark et al 



2004)) indicates a spatially flat universe with low 
matter density (around one-third of the critical 
density), the rest of the closure density is believed 
to be in an unknown form generically termed as 
the dark energy. It is this component that drives 
the late time acceleration of the expansion of the 
universe. 

The simplest possibility, which fits the data 
well, is that the acceleration is caused by the pres- 
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cosmology: large-scale structure of universe 



ence of a cosmological constant term, called the 
A term. When combined with the usual mat- 
ter term, the resultant model is referred to as 
the ACDM model. In this model the Hubble pa- 
rameter asymptotically approaches a constant at 
late times, thus causing the universe to acceler- 
ate. There are compelling theoretical reasons to 
believe that the dark energy density may not be 
a strict constant. There are several physical mod- 
els where dark energy is generated dynamically 
from an evolving scalar field (called quintessence) 
or even from alternate th eories of gravity (see 
Sahni k. Starobinskvl[200l for a recent review of 



models and methods of reconstruction of cosmic 
history) . Due to its simplicity, it has become pop- 
ular to phenomenologically model dark energy as 
an ideal fluid with an equation of state given by 
p = wp, where w is allowed to be negative. In this 
model w = —1 gives the usual cosmological con- 
stant. Models where w is a constant or a simple 
function of redshift have also been considered. 

Cosmological data is rapidly approaching a 
quality where we can start discriminating com- 
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peting models of dark energy. The effect of tiny 
departures from a strict cosmological constant on 
observations is sufficiently small to render such 
analysis unreliable if the nature of statistical noise 
in data is not well understood. Since the most de- 
tailed constraints on dark energy are derived from 
the luminosity distance to distant supernovae, 
we would like to be certain that their statistics 
i s well understood. Th e central limit theorem 
(|Kendall fc StuartJ Il977h ensures that, to a very 
good approximation, statistical noise due to a 
large number of random influences can be treated 
as Gaussian. There are several possible sources of 
no n-gaussianity in t he sup ernova data; according 
to iKolatt fc Lahavl (|200ll ) these could be 1) sta- 
tistical scatter due to location in the host galaxy 
and galaxy type 2) scatter due to dust absorption 
in the host galaxy, inter galactic medium or in our 
own galaxy or 3) due to lensing along the line of 
sight. Some of these processes are corrected for 
in the data reduction. It is also possible that the 
observed anisotropy is a result of collation of data 
of disparate quality, perhaps due to differences in 
seeing condition or in the data reduction process. 

Modern cosmolog y is based on the Cosmolog- 
ical Principle (CP) (jPeeblesI fl993h , which states 
that on the large scales the universe is statistically 
homogenous and isotropic. Even if supernovae 
were perfect indicator of distance, and if statis- 
tical noise in the supernova distances is Gaussian, 
an anisotropic universe could contribute a system- 
atic, direction dependent modulation in the data. 
Another possibility is that our galaxy could con- 
tain anisotropic, gray dust patches, and since the 
dust correction depends on reddening in the spec- 
trum, this sort of dust could remain uncorrected. 
Above arguments suggest that there is a strong 
case for investigating direction dependence in the 
supernova data. 

IKolatt fc Lahav (2001) have investigated the 
possibility of detecting cosmic anis otropy with 



79 hig h-z supernovae obtained fr om iRiess et al 



(jl998l ) and lPerlmutter et alj (| 19991) . In this paper 
we use a different statistic than used by them. We 
use the gold data set containing 157 supernovae 
(jRiess et al.1 120041 ) for our analysis. Although our 
methodology is derived from testing for isotropy, 
we shall use the term non-gaussianity interchange- 
ably. Our results can be interpreted as a system- 
atic directional dependence in the data due to any 



of the above mentioned possibilities. 

The outline of this paper is as follows. In |}2]we 
describe our methodology in detail. We present 
our results in JJ] and our conclusions in 2J 

2. Methodology: Extreme Value Statistic 

Throughout our analysis we have assumed a 
flat FRW universe. Since ACDM model fits 
the data quite well, we first o btain the best fit 
model to the full gold data set (jRiess et al.1 120041 ) 
and calculate the dispersion normalized residu- 
als Xi = [Mi ~ Ai ACDM (^; Om)]/^^), where the 
distance modulus fi = 5 log(rfi,/Mpc) -I- 25, the 
observed values being for a supernova at red- 
shift Zi, and cr M (zi) is the observed standard error. 
We shall consider subsets of the full data set to 
construct our statistic. We define the reduced x 2 
in terms of Xi as follows 
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where it should be noted that by 'reduced' we do 
not mean 'per degree of freedom', since we do not 
fit the model separately to the subsets of the full 
data. Here x 2 is an indicator of the statistical 
scatter of the subset away from the best fit ACDM 
model. 

If CP holds then the apparent magnitude of a 
supernova should not depend upon the direction 
in which it is observed but only on the cosmol- 
ogy. We divide the data into two hemispheres 
labeled by the direction vector rij, and take the 
difference of the reduced x 2 computed for the 
two hemispheres separately to obtain A% 2 . = 
Xnorth — Xsouth' wnere we have defined 'north' as 
that hemisphere towards which the direction vec- 
tor Hi points. 

We take the absolute value of Ax 2 . , since that 
is the quantity which determines the plane across 
which data is anisotropic, and then vary the di- 
rection n across the sky to obtain the maximum 
absolute difference 



A = 



cilAxlJ}- 



(2) 



To interpret our results we need to know what val- 
ues of A one might expect. The distribution of su- 
pernovae is not uniform on the sky, therefore, the 
number of supernovae in the two hemispheres for a 
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given direction varies with the direction n. There- 
fore one might expect the probability distribu- 
tion function -P(A) to be extre mely complicated, 
howe ver, extreme value theory (jKendall fc Stuart 
19771 ) shows that the distribution is, in fact, a sim- 
ple, two parameter Gumbel distribution, charac- 
teristic of extreme value distribution type I: 
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where the position parameter m and the scale pa- 
rameter s completely determine the distribution. 
To quantify departures from isotropy we need 
to know the theoretical distribution Pthcory(A). 
Since it is difficult to obtain it analytically, we 
have calculated it numerically by simulating sev- 
eral sets of Gaussian distributed Xi on the gold set 
supernova positions and obtaining A from each re- 
alization. We plot this distribution in Fig [1] as the 
broken curve. We find that the distribution closely 
resembles a Gumbel distribution. 

If the data do have directional anisotropy then 
an independent possible test for non-gaussianity 
is obtained by constructing the bootstrap distri- 
bution Pbs(A) in the following manner. The ob- 
served Xi are assumed to be drawn from some 
unknown, direction dependent probability distri- 
bution. We shuffle the data values Zj, m(zi) 
and a m (zi) over the supernovae positions, thus 
destroying any directional alignment they might 
have had due to anisotropy. We show in the next 
section that this gives us yet another way of quan- 
tifying non-gaussianity. 

3. Results 

Our main result is plotted in the Fig [1] We 
find that the theoretical distribution P t h e ory(A) 
(broken line) assuming Gaussian distributed Xi 
indicates that our universe has about one sigma 
smaller anisotropy than the mean of the distri- 
bution. However, as mentioned in the last sec- 
tion, if the residuals fa are non-Gaussian then a 
more appropriate estimate of departure from non- 
gaussianity would be the bootstrap distribution 
P BS (A). We have plotted P B s(A) obtained by 
randomly shuffling Xi on the given supernovae po- 
sitions in the same figure. We find that the ob- 
served value of A is more than one sigma away 
from the mean of this distribution as well. 



One problem with the bootstrap distribution is 
that we expect it to be shifted slightly to the left 
of the theoretical distribution. This is due to the 
fact that theoretical distribution is obtained by as- 
suming Xi t° be Gaussian random variates with a 
zero mean and unit variance. Therefore theoret- 
ical Xi s are unbounded. However, the bootstrap 
distribution is obtained by shuffling through a spe- 
cific realization of Xi > an d they have a maximum 
value such that \xi\ < Xo- Since the bootstrap re- 
alizations have bounded Xij they should produce, 
on the average, slightly smaller values of A as 
compared to what one expects from a Gaussian 
distributed Xi ■ For a large number of supernovae 
this bias is expected to vanish. Since the gold data 
set contains only 157 supernovae, another concern 
is that our statistics may not be reliable enough 
due to a lack of uniform sky coverage. We have 
made a few checks to investigate these concerns: 
We simulate the same number of supernovae as in 
the gold data set with randomly chosen positions 
on the sky. We then process the simulated data in 
exactly the same way as the actual data. A typical 
result is shown in FigO There are a few things to 
be noted. 1) The simulated data gives a A that is 
consistent with the theoretical distribution, indi- 
cating that a uniform sky coverage is not a strict 
requirement in this statistic. 2) The theoretical 
and bootstrap distributions for the simulated data 
look identical in shape and 3) the peak of the boot- 
strap distribution is shifted leftwards, as discussed 
above. This shift is on the order of 10 per cent. In 
Fig Q] we show that the bootstrap distribution has 
a mean that is about 40 per cent shifted from the 
theoretical distribution, and has a different shape, 
independently indicating non-gaussianity. The ex- 
cess shift cannot be reconciled with the theoretical 
distribution by a simple scaling of the error bars. 
To produce a rightward shift we would need to 
increase A, which can be done by decreasing the 
error bars on supernovae by a constant scale fac- 
tor. However, this would also produce a larger xi 
for the best fit ACDM model. The data actually 
gives xl — 1-14, and scaling would make it larger, 
thus making the primary fit worse. 

We find that the gold set is maximally asym- 
metric in the direction (I = 100°, b — 45°). We 
designate the two hemispheres as 'hot' or 'cold' ac- 
cording to largeness or smallness of their reduced 
X 2 (as given in Eq. [IJ with respect to the best fit 
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ACDM model. In Fig [3] we have plotted the resid- 
uals Ayit = /Xdata — macdm for the hot and cold 
subsets of supernovae. The zero line is the best fit 
ACDM model. As expected, the cold supernovae 
show slightly smaller dispersion compared to the 
hot ones. Figure [3] also shows that the ACDM 
model seems to fit the hot and cold supernovae 
equally well. This is borne out by a parameter 
estimation, which is tabulated in Table [TJ where 
we find that the best fit ACDM model for hot su- 
pernovae gives fi m = 0.30 and the cold ones give 
f2 ra = 0.31, so the difference is only a few per 
cent. However, the situation is not the same for 
model where we have assumed a constant equa- 
tion of state p = wp. We find that the model 
parameters for the hot supernovae in this model 
are f2 m = 0.51 with w — —4.53 and for the cold 
supernovae f2 m = 0.32 and w = —1.03. The value 
of the Hubble constant is relatively quite robust, 
showing that most of the effect is coming from 
high-z supernovae. The large difference in the val- 
ues for the constant w model shows that the level 
of non-gaussianity indicates that constraints on a 
more complicated dark energy model are not as 
ro bust. Perhaps this explains the intriguing result 



Alam et al. ( 20061 ). that the data seems to fit a 



ACDM model as well as a model with a strongly 
evolving dark energy. 

4. Conclusions 

We have used the extreme value statistics on 
the gold data set and shown the presence of non- 
Gaussian features in the data. We find that there 
is a direction of maximal asymmetry in the data 
across which data seems to imply different cos- 
mological models for a constant equation of state, 
although, the constraints on the ACDM model 
are found to be robust. We have discussed how 
this could either imply non-gaussianity in the data 
due to various possible physical processes or due 
to anisotropy in the universe. Our results can- 
not be trivially understood by scaling the error 
bars. Since this analysis has been done within 
the framework of an FRW model, it is obviously 
very difficult to quantify the precise meaning of 
this anisotropy. We have discussed that we need 
to be very careful in interpreting dark energy be- 
yond the cosmological constant model since it is 
possible that systematic noise may masquerade as 
evolving dark energy. Further work is required 



to fully understand the statistic that we have in- 
troduced in our analysis and will be discussed in 
greater detail in a future publication. 

Shashikant thanks CSIR for providing financial 
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Fig. 1. — The solid line shows the bootstrap prob- 
ability distribution Pes (A) obtained by shuffling 
supernovae data on the sky while keeping the su- 
pernovae positions fixed. The broken line shows 
the theoretical distribution P t hcory(A) obtained by 
simulating Gaussian distributed x%, as described 
in the text. The vertical lines show one and two 
sigma regions. The observed value for A for the 
gold data set is seen to be around one sigma away 
from either distribution. 



Table 1: The model parameters for the hot and 
cold set are tabulated here. WCDM refers to a 
model where the dark energy has a constant equa- 
tion of state. 
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Fig. 2. — We plot the result of a typical simu- 
lation where 157 supernovae positions were gen- 
erated randomly and populated with Gaussian 
noise. Similar to Fig [1] this plot shows the Gaus- 
sian vs. the bootstrap distribution. The simulated 
universe is seen to be consistent with the theoret- 
ical distribution, thus indicating that our statistic 
does not require a uniform sky coverage 
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Fig. 3. — The top panel shows the distribution of 
cool supernovae and the bottom one shows the hot 
supernovae. Horizontal line is the ACDM model 
in both panels. The total number of supernovae 
in the top panel is 74 and the bottom panel 83. 
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